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Abstract 

The general local, nondissipative equations of motion for a quantized vortex 
moving in an uncharged laboratory superfluid are derived from a relativistic, 
co-ordinate invariant framework, having vortices as its elementary objects in 
the form of stable topological excitations. This derivation is carried out for a 
pure superfluid with isotropic gap at the absolute zero of temperature, on the 
level of a hydrodynamic, collective co-ordinate description. In the formalism, 
we use as fundamental ingredients that particle number as well as vorticity 
are conserved, and that the fluid is perfect. No assumptions are involved 
as regards the dynamical behaviour of the order parameter. The interaction 
of the vortex with the background fluid, representing the Magnus force, and 
with itself via phonons, giving rise to the hydrodynamic vortex mass, are 
separated. For a description of the motion of the vortex in a dense laboratory 
superfluid like helium II, two limits have to be considered: The nonrelativistic 
limit for the superfluid background is taken, and the motion of the vortex is 
restricted to velocities much less than the speed of sound. The canonical 
structure of vortex motion in terms of the collective co-ordinate is used for 
the quantization of this motion. 
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I. INTRODUCTION 



Vortices are the fundamental line-like excitations of a fluid. Their motion governs a large 
number of properties of the fluid under given external conditions. Superffuids are peculiar 
in that the vortices can exist only in certain classes, distinguished by integers. This fact 
that vorticity can only arise quantized, results in far-reaching consequences for the character 
of vortex generation and motion in quantum fluids, as compared to simple perfect fluids. 
The aim of the present investigation lies in the derivation of the general local equations of 
motion for a quantized vortex in a laboratory superfluid at absolute zero, thereby starting 
from a minimum of very basic assumptions about the co-ordinate invariant nature of the 
hydrodynamic conservation laws in a relativistic perfect fluid and the Magnus force acting 
on the vortex. In particular, no assumptions will be made with respect to the underlying 
dynamics of the order parameter of the superfluid. 

The most fundamental property of a vortex living in a superfluid is its topological sta- 
bility. It is a topological line defect belonging to the first homotopy group Ui(G/H), which 
contains the equivalence classes of loops around a singularity: The circle S 1 in real space is 
mapped onto the coset space G/H, and the equivalence classes identified. Here, G is the 
symmetry group of the system under consideration, and H is the little group, a subgroup 
of G, characterizing the partial symmetry of the system which remains after the symmetry 
of the larger group G has been broken. This coset space represents the manifold of the 
degenerate vacuum states after spontaneous symmetry breakdown, and can be identified 
with the order parameter. In this paper, we will deal exclusively with a very simple such 
coset space, namely the circle S 1 (the group G is U(l) and H is the identity), so that the 
mapping in question is between the circle in ordinary and that in order parameter space. 
The winding number N v G Z = II 1 (S' 1 ) classifying the paths around the line equivalent to 
each other, indicates the number of times one is 'winding' around the singular line in order 
parameter space, if one is going once around it in real space (see Fig. [l]). The complex order 
parameter of a spontaneously broken U(l) symmetry, with phase 6, then gives rise to the 
quantization of the nonrelativistic superfluid circulation T s in units of k — h/m, 

T s = I v s . ds = — lv9 ■ ds = N v - — = N v - k, (1) 
J m J m 

provided the superfluid velocity is identified as v s = (h/m)V9. For a vortex in the superfluid 
helium II, the superfluid particle's mass m = 6.64 • 10 _27 kg, so that k ~ 10~ 7 m 2 /s. The 
definition in ([I]) is the classical kinematical definition of circulation We will see below 
that a relativistically invariant, dynamical definition has to use the line integral of the 
momentum, rather than that of the velocity, to gain truly invariant meaning. However, 
in any case, be it relativistic or nonrelativistic, the stability (and very definability) of the 
quantized vortex as a fundamental object able to persist while moving in the superfluid stems 
essentially from its topological properties. The superfluid cannot be continously deformed 
into a vortex-free state. On the other hand, this feature of topological stability has its limits 
if large quantum statistical fluctuations of the order parameter become important. In the 
dense, unpaired superfluid helium II, they become large on the same coherence length scales, 
on which the order parameter modulus varies, of the order of the inter-particle distance. 
Hence a vortex, having curvature scales of the singular line, which are of order the inter- 
particle distance, is not definable topologically and thus does not exist as a well-defined line 
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defect in order parameter space. On these scales, it has to be defined using the full quantum 
many-body structure of the dense superfluid and the corresponding vortex wave functions. 

This hydrodynamic level of accuracy can only be afforded by the fully relativistic de- 
scription of the dense liquid we present in the section which follows, too, because the same 
problem of incorporating microscopic many-body structure excitations into the description 
remains, which exists already in the nonrelativistic superfluid. A hydrodynamic formalism 
is thus a necessity to describe vortex motion properly, as long as this problem of microscopic 
dynamics has not been solved. The covariant treatment, in turn, provides structural insights 
which are not given by more conventional treatments. 
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FIG. 1. A rectilinear quantized vortex with negative winding number N v = —1. The superflow 
in the x-y plane, indicated by the arrows, is in clockwise direction along the gradient of 6. For a 
vortex with positive winding number the superflow is anticlockwise, correspondingly. The branch 
cut is chosen to be at = ±7r and the circulation vector is pointing into the figure. 



II. DUALITY OF VORTICES AND CHARGED STRINGS 

We establish in this section the properties of a vortex as a stringlike fundamental object. 
For this purpose, the means of dual transformation will be instrumental (where dual is 
meant in the sense of being 'equivalent as a physical system'). We will introduce the duality 
of interest here in its relativistic context (established in Refs. [|3]|4|], elaborations in various 
directions can be found in Refs. 0- fl2fl), subsequently reducing it to its nonrelativistic 
limit. 

Conventions are that the velocity of light c is set equal to unity and the signature 
r}fj, v = diag (—1, 1, 1, 1) of the metric is being employed. Greek indices mean spacetime indices 
and take the values 0, 1, 2, 3. The sign convention we use for the Levi-Civita pseudotensor 
is e 0123 = +1 = —6oi23 in a Lorentz frame. The dual of a p-form / in four- dimensional 
spacetime has the contravariant components yPi-^-p = (l/p\)f ai ,^ ap e ai '" apl3l '" f34 ~ p . 

The hydrodynamic conservation law we will use as a first basic ingredient of our theory for 
an uncharged relativistic fluid is that of particle number. The most familiar mathematical 
form of nonrelativistic number conservation is provided by d t p + div(piT s ) = 0. This is 
expressed covariantly, using p— forms [ITT 



as 



dA*j = (d = dx a d a ). (2) 

The number current one-form is defined j = pu = pu a dx a , where dx^ is a one-form basis 
dual to a co-ordinate basis and p the rest frame number density. We omit the subscript s 
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from the relativistic quantity u, whose normalisation will be required by its interpretation 
as a four- velocity to be u a u a — — 1. 

The momentum of a perfect fluid particle is defined as the one-form p, with components 

El 

Pa = I^U a . (3) 

The chemical potential p = de/dp, where e is the rest frame energy density, then plays 
the role of an effective rest mass in a Hamiltonian quadratic in the four-momentum. The 
momentum and number current density one-forms are thus related through 

£.-—-' (A\ 
Pa ,_o Pa Ja ■ V*) 

The quantity K = h 2 (p/ p) is the stiffness coefficient against variations of the order parameter 
phase (a quantity to be defined below) in the free energy, cf. [ll| and section 3.1 in [20 



The speed of sound c s is related to the quantities introduced in the above, for a barotropic 
fluid @, as c 2 = d(\np)/d(\np) = (K/h 2 )d 2 e/dp 2 . 



We define the vorticity as the dual of the exterior derivative of p: 

uj = *d A p <^> *uj = -d A p . (5) 

The vorticity thus defined is required to be conserved, in the form 

dA*u = O dAdAp = 0. (6) 

Outside the cores of the vortices, where the quantity vorticity uj defined above equals 
zero, the momentum of the fluid always remains proportional to the (exterior) derivative of 
a scalar 6, which we identify with the phase of the U(l) order parameter pertaining to the 
superfluid: 

Pa = hdJ. (7) 

The circulation is in the relativistic case defined as the integral of the momentum. This 
choice of definition stems from the fact that a particle's rest mass, occuring in the nonrela- 
tivistic definition of the circulation ([[]), is an undefined quantity in a fully relativistic, dense 
superfluid. The proper circulation thus reads 

% = j P^dx** = N v h , (8) 

and is quantized into multiples of Planck's quantum of action: The constant of propor- 
tionality h = h/2ir in (|7|) stems from the covariant Bohr- Sommerf eld quantization of 7 S , 
i.e. from the fact that we require p to be a quantum variable. It may be worthwhile to 
point out in this context that particularly striking evidence for the fundament ality of the 
one- form quantity momentum, p = p^dx^, is given by consideration of superfluid rotation 
in the presence of gravity [15], e.g., in neutron stars. Superfluid irrotationality is equivalent, 



for axial symmetry, to p^ = 0, whereas the contravariant axial component pr and thus the 
contravariant axial velocity do not have to vanish. 
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The fact that quantization of circulation as a line integral of momentum is crucial for 
a proper understanding of the nature of the line defect vortex is not connected to the fact 
that the superffuid is relativistic or nonrelativistic. The point is rather that the vortex is 
in general nothing but a singular line of zeroes in the order parameter manifold, designated 
by a single number N v , which represents the order of the pole. This fundamental, intrinsic 
property of the vortex is not related to a quantity like mass, characterizing the matter with 
which it interacts, but is independent of matter properties. 

The normalisation of u^, by using the metric and writing u a u a = —1, is promoting 
into a four- velocity. It leads, with (|3|) and ([7p, to the relation 

p = h{-d^d^e) 1/2 (9) 

between the phase and the chemical potential. This is a relativistically covariant version of 
the Josephson equation faQ = — (/i + ^mv 2 ), familiar from nonrelativistic condensed matter 
physics |^3||. The Josephson equation expresses, generally, the conjugateness of the two 
canonical variables phase and number density in the hydrodynamic limit. 

The dual of the current j = j^dx^ is a 3-form *j = j^e^^ dx u A dx a A dx^, which is, 
according to the conservation law above, closed. We define the field strength H by 

*j = H. (10) 

The field H^ ua is totally antisymmetric in its three indices and has, by definition, only four 
independent components. In a simply connected region, H is exact, i.e., it is the exterior 
derivative of a gauge 2-form b = b^ v dx^ A dx u : 

Ha/3-y = dabp-y + d^ba/3 + dpb^ a , (11) 

The field strength H is invariant under gauge transformations b — > b + d A A, where 
A = A a dx a is an arbitrary 1-form. Thus 

b>xv — ► fyui/ + d^Ky — dyk^ =^ Hixvcl —> H^va . (12) 

We then have the following sequence of relations, which define the correspondences of p, p, 6 
and b: 

*j = —*p = h— *d6 = d Ab = H, (13) 
p p 

The dual transformation in the conventional sense obtains if we neglect any possible 
variations of p/p = K/h 2 in the spacetime domain of interest. Then, b has only one degree 
of freedom corresponding to the order parameter phase 9. 
From the equation H = (p/p)*p, it follows that 

*dA (j* H ^ = * d A = *dAp = u, (14) 

which is the general field equation for the motion of vortices with conserved vorticity u). 
It is a field equation analogous to the inhomogeneous Maxwell equation. In case that the 
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stiffness ratio of density and chemical potential is a constant Kq in space and time, and in 
Lorenz gauge *d A*b = 0, the wave equation of the gauge field takes the familiar form 

Db=-(K /h 2 )u J , (15) 

where the d'Alembertian □ = . The homogeneous Maxwell equation for strings is, by 
definition, for any value of p/p, 

d A H = (arbitrary p/p = K/h 2 ) . (16) 

The advantage of using b over 8 consists in the fact that it can be chosen to be single- 
valued, which is to be contrasted with the necessary multivaluedness of 6 in the presence 
of a vortex. The physical significance of b is that it represents a generalization of the 
stream function concept encountered in classical nonrelativistic hydrodynamics [0. Namely, 
integrating 

/ h^dx* 1 A dx v = / pu x ex^dx^ A dx v A dx a , (17) 

we see that the integral of b over any 2-surface enclosing a 3-surface in spacetime 
is given by the superflow flux through the 3-surface (the symbol | ■ • • | indicates ordering of 
the indices contained between the vertical lines with increasing numerical value). The 2- 
form quantity b can thus be called a generalized stream tensor. The relation flT7j) expresses 
a generalization of the usual concept of stream function (for which the incompressibility 
condition p = po holds), to flows for which the density can vary. 

In particular, in the nonrelativistic case, which has u° = 1, setting p = po, 

- I bji dx i A dx j = p V( 3)n , (18) 

that is, the specified surface integral of the purely spatial part of b is the spatial volume 
enclosed by the surface ^dVL times the density po, and thus the bulk number of particles 
contained in the volume ^fl. 

The superfluid we will be dealing with later in section |T| is nonrelativistic. We thus have 
to rewrite (|13D in a form having only Galilean invariance. This is accomplished by writing 
(O) (with p = po) for spatial and temporal indices separately, 

PoU°eoijk = dkhj + djbki + dibjk = Hij k 

PoU l € i0jk = d k b j — djb ok + -d t b jk = H Qjk . (19) 

and taking the Galilean limit c — > oo (u° — > c, u l — >• v\ ; we have temporarily reinstated 

the velocity of light for this purpose). We then get the relations (in which e^jk = —tijk = 
-\fgn ijk )'- 

- y/g n ijk p = d k bij + djb ki + dibjk , (20) 

\/g nijk v l s = d k ipj - djipk , (21) 
where the vectorial version of the usual stream function M is defined through 

A = WPo • (22) 

The determinant of the spatial co-ordinate system we are using is designated g, and n^k = 
nujk] = ±1 is the unit antisymmetric symbol. 
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A. The Magnus force 



The fundamental force acting on a superfTuid vortex at the absolute zero of temperature 
is the Magnus force. This force is now shown to be equivalent to a stringy generalization of 
the Lorentz force. Given this identification, we will be able to write down electrodynamic 
correspondences in the next subsection. 

The Magnus force acting on a vortex in a nonrelativistic superfluid has the standard 
form 

F M = mp f s x \X - v s j 

= lsPo X' x (X - tf 5 ) , (23) 

where v s is the superflow velocity far from the vortex center located at X l (t, a) (i.e. the flow 
field at the vortex location without the contribution of the vortex itself), and a is the arc 
length parameter labeling points on the vortex string. The circulation vector T s = N v kX' 
of the vortex is, for positive Ny, pointing along the z-direction in a right-handed system. 

In the second line, the force is written in a form manifestly independent of a 'mass' value. 
A relativistic generalization of the above expression then reads 

(F M ) a = ls H aiiU X»X w . (24) 

We introduced the relativistic string co-ordinates X^ = X M (r, cr) , and let X ,fl = dX^/da 
be the line tangent, as well as X^ = dX^/dr the vortex velocity. The essential features of 
the vortex as a two-dimensional object, living in spacetime, are represented in Figure 0. 




FIG. 2. The vortex is represented by the string world sheet embedded in spacetime and hence 
described by the co-ordinates X^(r, a). The tangent space basis vectors have components in 
the timelike and X'^ in the spacelike direction on the world sheet, where the metric is given by 

lab = g»v dx^/dc dx»/d( b , with a, b = 1, 2, c 1 = r, C 2 = <?■ 



s 



It is to be stressed that both of the forces (|23|) and its relativistic counterpart (|24"D are 
forces per unit cx-length and are of topological origin. They do not depend on the local 
shape of a line segment, represented in the relativistic case by the world sheet metric 7 a b, 
but only on the local external field H afJjU , generated by other line segments and an externally 
imposed flow field (cf. the equation for the corresponding action fl25|) and the discussion at 



the beginning of section I I 1 B 



The Magnus force can be derived from the variation of the action (cf., in particular, the 
Refs. 011,0) 

S M = ls Jdrdab^X w 

= 7s ldSb^\ (25) 



that is to say 



5S 



M 



5X a 

The singular vorticity tensor components are thence given by 



ls H ailv X»X w = (F M ) a . (26) 



^ = Is 



J j drda (x*X' v - X V X'^ 5 {A \x - X(r, a)) . (27) 



From the vorticity tensor, we obtain the components of the usual vorticity vector by choosing 
X° = t, a choice possible in any global Lorentz frame fTT|: 



u 0i = ls J daX'^W (x - X(t, a)) . (28) 

For a rectilinear line in ^-direction, uj 0z = 7 s 5(x — X)5(y — Y), which integrated over the 
x-y plane gives the circulation. The vorticity tensor components (^) are, then, to be 
understood as a generalization of a quantity surface density of circulation, in the sense that 
ls = §P = JJdAp = -ff*u> = JJ -e lfluM uj^dx a A dx? = JJ u^d 2 S\^\. Integration 
is over a (sufficiently small, i.e. local) 2-surface with element d 2 S^ v = —^^ v \ a p\ dx a A dx 13 , 
which is threaded by the vortex world sheet. 



B. Electromagnetism of vortex strings 

It proves very suggestive to cast the laws we found into the language of Maxwell's elec- 
trodynamics. From the familiar Lorentz force law 

(-F'Lorciitz)a = ^F^X^ , (29) 

one can define that the 'electric' and 'magnetic' fields are represented by a projection of 
on the local vortex axis, i.e. on the tangent space vector along the spacelike direction 



on the world sheet 141: 
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= H^pX' 13 — e^ ua pj a X'^ , 
E = Po uxX' , B = e 0123 p u°X'. (30) 

The second line is valid for a global Lorentz frame, in which, as already mentioned above, 
we can choose X° = t. In general, the density p in the field H^p, and thus also F^, is 
arbitrary; at the location of the line itself and for the Lorentz/Magnus force law, however, 
we set p = Pq. In our sign convention for €0123 = — 1, the local 'magnetic field' is antiparallel 
to the local tangent. 

The local 'electromagnetic' 3-potential will accordingly be defined via 

a^ = b^ = b^ v X lv . (31) 

In an arrangement of cylindrical symmetry, fulfilled by a ring vortex, the 3-potential 
is the vector a M = (a°,a r ,a z ) = {—bQ^b r ^b z ^). The component ips = &o</>/Po corresponds 
for stationary flows and p = po to Stokes' stream function 0j. The function ips fulfills 
27ii/js = 27r J r(u r (r, z)dz — u z (r, z)dr), according to (pT7f) . The static scalar potential of our 
'vortodynamics' is thus in the cyhndrically symmetric case given by a = po J r(u z dr — 
u r dz). Integrating this relation to obtain the value of a°(r ,z ) at some point r Q ,zo in the 
superfluid, the line integrals run over the envelope of an arbitrary surface of revolution, 
which is generated by rotating the integration line, joining the point at r , Zq and a point 
on the z-axis of symmetry ||. 

In making the above definitions (|3"0|), we identified q = 7 S = N v h with the 'charge' 
per unit length of vortex. This vortex 'charge' is thus of topological nature, quantized by 
the index of the homotopy group member (in units with h = 1, the charge q = 2irN v ). 
The quantization of the charge vanishes in the purely classical limit of h — > 0, whereas the 
'electromagnetic' field strength F is understood to be derived from a (superfluid or not) 
conserved number current. 

We can define the vortex current four-vector in a fashion analogous to the definition of 
the potential in (|3TD, from the minimal coupling term represented in (p5|): 

W = = cu^Xl , (32) 

that is, we project the vorticity current on the local vortex axis. This completes the picture 
of the electromagnetism analogy we wanted to develop. We have constructed, within the 
semiclassical realm, a fundamental topological object, the vortex, which moves according to 
familiar laws of electrodynamics, by projecting the string equations of motion on the local 
spacelike vortex tangent axis. The field equations have the form of Maxwell equations local 
in a, 

d A F = , *d A (h 2 *F/K) = O , (33) 

where the field F is the usual 2-form and the current O a 1-form, whose components are 
defined in ([ft]) and (|3"2"|), respectively. The vortices can be understood as fundamental, 
'charged' objects, generating 'electromagnetic' fields, which in turn act on them by the local 
Lorentz force in (f29]).[] 



1 We remark that we have chosen not to incorporate the factor 4-7T, appearing in the cgs units 
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The equations (|33| ) are representations of the conservation laws of particle number, vor- 
ticity, and the equation (||) characterizing properties of the medium, which relates these 
conservation equations to each other. The homogeneous 'Maxwell' equation is a projection 
of number conservation, whereas the inhomogeneous 'Maxwell' equation is equivalent to a 
projection of (the conservation of) vorticity, combined with a property of the medium. Vice 
versa, we conclude that the Maxwell equations of ordinary electromagnetism can be cast 
into the form of conservation equations of relativistic perfect fluid hydrodynamics, provided 
we admit the identifications specified in ([13]). 



C. Comparison of relativistic and nonrelativistic notions 



The identifications we have made pertain to the vortex as a stringlike fundamental object 
and are suitable for a relativistic framework. In the remainder of this paper, however, 
we will only be concerned with nonrelativistic Eulerian flows.0 At this stage it is thus 
useful to compare and contrast the relativistic notions and those of classical nonrelativistic 
hydrodynamics 

The most important difference of the two frameworks occurs if we consider the two 
definitions (|T|) and (^) for the circulation. The kinematical quantity velocity suitable for 
its definition in nonrelativistic circumstances can not be used for its relativistic definition. 
The dynamical quantity average angular momentum per particle, though, which is quantized 
into units of h because of the existence of a macroscopic superfluid phase, is a well-defined 
quantity for any superfluid. Because of this fact, the definition in @ is the truly invariant 
definition of circulation. The only possibility to use a 'velocity line integral' for the definition 
of circulation is to employ the perfect fluid relation p a = fiu a . Then, however, one has the 
complication of a possible position dependence of /i. 

A prescription to translate from the fully relativistic superfluid to terms of standard 
nonrelativistic hydrodynamics is afforded by the limits fi — > m, u a — > (l,v % s ), and the 
replacement 7 S — > T s . 



III. CANONICAL STRUCTURE OF VORTEX MOTION 



A. Phonons and the Vortex Mass 



In quantum electrodynamics, the fundamental entities interact between each other and 
with themselves by photons. The corresponding quanta in our view of a superfluid are 



Maxwell equations, into the definition of the vortex charge, i.e. wrote the 'Maxwell' equations in 



Heaviside-Lorentz units |22|. 



2 The term 'Eulerian' is used for a fluid obeying in its (non-dissipative) dynamics a continuum 
version of Newton's equation of motion, the Euler equation (there exist also generalizations of this 
equation for curved space-time backgrounds, cf. §22.3 in p" 
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phonons, the excitations of the 'vacuum ether' surrounding the 'charged' strings. The role 
of these excitations is addressed in this section. 

The vortices acquire a nonzero hydrodynamic mass from their self-interaction with these 
phonon excitations. This is most easily understood if one realizes that in a 2+ld superfluid 
there exists a direct correspondence between the genuinely relativistic electron-positron pair 
interaction via photons in quantum electrodynamics and a 'relativistic' point-vortex-point- 
antivortex pair interaction via phonons P4] , |2"5| (we will use throughout inverted commas to 
distinguish pseudorelativistic = 'relativistic' behaviour from the actual Lorentz invariance 
of section |TJ). In the light of the vortex — fundamental object correspondence expounded in 
the preceding section this is only natural: The only 'vacuum' excitations ('vacuum' in the 
sense of the spontaneously symmetry-broken superfluid vacuum) , by means of which different 
vortex line segments can interact are, in the hydrodynamic limit, phonons. What remained 
to be done is to convert the quite direct 2+ld correspondence into a 3+ Id correspondence 
by defining the respective quantities local on the string, i.e. as functions of a. Then, the 
phonons give 'relativistic' fields propagating on a nonrelativistic background of Eulerian 
superflow, and mediating interactions between (singular) vortex line segments. 

To explain this further, we begin by considering that the unit circulation vortex carries, 
per unit length, the hydrodynamic self-energy 



F h2p0 
-G/self — -; 

Arrm 



, , 8R, 
In 



c 



(34) 



with it. The energy E se n is the energy of a vortex sitting at some fixed place, that is, its 
rest energy in its rest frame. 

The infrared cutoff in the logarithm is in the static limit equal to the mean distance of 
line elements, respectively, in the localized self induction approximation, proportional to the 



local curvature radius R c of the line. The constant C parameterizes the core structure |26 
and has order unity. In classical hydrodynamics, with a hollow core of radius £, C — 2, 
whereas in a model of constant core vorticity one has C = 7/4. In the Gross-Pitaevskii 
framework ( |27|], p5[|), for singular vorticity, it turns out that C = 1.615 p6 |. 
If we separate off in the order parameter phase two parts, 

9 = @bg + 9 P h , (35) 

a part due to a background flow 9b g , and a part which describes sound excitations 9 ph <C 9b g 
on this background, we can carry through the program of the last sections for 9 bg and 9 p h 
separately. For the phase 9b g , we can not retain the Lorentz invariance of the equations we 
derived there (our fluid is very much Eulerian), and have to use the Galilei invariant set of 
equations fl2"0|)-(]2~TD. On the other hand, for the part 9 p h, a 'relativistic' equation from which 
to start the dual transformation derivation for 9 P h is the wave equation 



d^e ph 



2 



- — 



9 ph = . (36) 



This is 'relativistic' in the sense that the speed of light c is replaced by the speed of sound c s 
in a Lorentz invariant scalar wave equation for 6 p h. From this part of the phase we can derive 
'electromagnetic' phonon wave field strengths obeying the 'relativistic' Maxwell equations 
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The nonlinear equation of motion of the superfTuid, the Euler equation, gives in its 
linearized version, together with the continuity equation, the above equation of motion (|36D 
for # ph -C 9b g , if the background flow is at a velocity much less than that of sound. For general 
background flows with velocities {h/m)\ V6 bg \ < c s , there results a scalar wave equation for 
sound in a curved Lorentzian signature background metric with nonzero curvature. See the 
nice discussion in |3l| (also cf. section IV in |HJ). We consider, however, 'nonrelativistic' 



background flows of small Mach number, so that these corrections are not of importance here 
(though they may be important in other contexts like for the Aharonov-Bohm interferences 
of phonons leading to the Iordanskh force J32[|), and the 'acoustic metric' is taken to be that 
of a global Lorentz frame. 

We now have defined the entities vortices as objects defined by the superfluid vacuum 
background and 'relativistically' interacting by the small density and phase perturbations 
in this vacuum medium. Hence, they obey the Einsteinian mass-energy relation 

M c 2 s =E seU . (37) 

The mass Mo is the hydrodynamic vortex mass in the rest frame of a vortex in equilibrium. 
If we consider a neutral, unpaired superfluid, this mass is dominant compared to other 
possible sources as the so-called backflow mass or the core mass corresponding to the normal 
fluid in the core This is true because we are in the limit of £ <C R c , in which the 

logarithmic divergence of the hydrodynamic mass dominates other possible contributions. In 
field theoretical terms, the energy contribution of the Goldstone boson field (i.e. the phonon 
field) is not screened by the gauge field like in superfluids with a local dynamical gauge field, 
which are charged in a conventional sense. This leads to the logarithmic divergence of the 
energy associated with the spontaneously broken global symmetry. 

It is useful to compare the definition of the vortex mass above to that of the classical 
electron radius r c in electromagnetism. For a homogeneously charged electron of radius r c , 
this definition reads m e c 2 = (3/5)e 2 /r c , where the factor 3/5 stems from the electron ('core') 
structure, which was assumed to be homogeneous. In the case of the electron, the mass serves 
to define, via the expression for the field energy, the value of the classical electron radius. 
In the case of the vortex (a line 'charge'), the field energy of flow around the vortex serves 
to define the renormalized, hydrodynamic vortex mass and a cutoff needs to be introduced 
because this field energy is logarithmically divergent. 

The knowledge of the mass M will be sufficient for our purposes, as the formalism 
laid down in these pages should apply for the actual, dense superfluid and consequently is 
restricted to velocities of the vortex and the background superfluid much less than the speed 
of sound, indeed even much less than the Landau critical velocity of roton creation. Also, we 
will consider scales much larger than £. Correspondingly, the (local) frequencies of vortex 
motion will be much less than the typical frequency uo s = c s /£. Scaling this frequency with 
parameters appropriate for helium II gives 

„. = 1 = 9.4. 10"* (38) 

which turns out be very close to the roton minimum (a> ro to n — 1.13- 10 12 Hz at p ~ 1 bar). We 
have scaled the coherence length with the Lennard Jones parameter of the helium interaction, 
&lj = 2.556 A, and assume is approximately unity. The condition of 'nonrelativistic' 
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velocities on scales much larger than £ is then equivalent to lu <C uj s . For 'relativistic' vortex 
motion frequencies in the order of and above u s , the vortex mass will be frequency and wave 
vector dependent pSfl. These high frequencies, employed in a semiclassical treatment, may 
be of interest in dilute superffuids. In what follows, we however restrict ourselves to describe 
'nonrelativistic' vortex motion, which is the case of interest in a dense superfluid. 

In the context of the question of vortex mass, it is advisable to point out that it is 
unnecessary, and indeed misleading, to use any postulated dynamical behaviour for the order 
parameter in the Lagrangian form (for example, of the Gross-Pitaevskii or time dependent 
Ginzburg-Landau variety). If we discuss the proper way of deriving the mass fl37p in a 
dense superfluid, we just need that the phase moves according to the Josephson relation 
m = -(fj, + \mvf) ||29|| . The validity of the Josephson relation, in turn, is not dependent 



on a Lagrangian for the order parameter, and is derivable from (Hamiltonian) superfluid 
hydrodynamics alone |30[|. 



B. The Vortex Lagrangian 

We have argued that the vortex has mass. Being a line-like, stable topological object, we 
assume it further to develop restoring forces if it is deformed, starting from some equilibrium 
position. That is, the vortex is assumed to have elastic energy, arising from the (local) 
interaction of a particular line segment with adjacent line segments. 

Let us first spend a few words on the relation of the self-action we will write down to 
'relativistic' vortex motion properties. In the 'nonrelativistic' case under study, the dif- 
ferential string arc length will be written as \f^(t, a)da (we remind the reader that the 
inverted commas are used to distinguish pseudorelativistic from proper relativistic terms). 
The self-action of the singular Nambu stringf] is proportional to the world sheet area p6l : 



^Nambu = -^T V~l(C°, C) A , (39) 



which is given by the double integral of the square root of the negative world sheet metric 
determinant y/— 7 over the timelike (° = r and spacelike C 1 = o co-ordinates, parameterizing 
the world sheet (cf. Figure ||]). The Nambu action is minus the constant string tension fiT 



(i.e., a constant self energy per unit length), multiplied by this area \2l\, 56|,P7|. The writing 
for the string arc length we will employ thus indicates that we are dealing with a 'Galilean 
limit' of the 'world sheet area' interval, in which the only metric element is the proper length 
interval of the line. 

The Nambu string is structureless (i.e., the core extension is negligibly small on the scales 
of interest), and has constant string tension. The actual superfluid vortex has a (quantum 



3 We refer to this elementary extended object for brevity simply as 'Nambu' string, though it has 
been introduced independently by Nambu and Goto |34 , [3q] , and thus is (properly) also quite often 
referred to as 'Nambu-Goto' string. 
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mechanical) core structure of extension £ which is comparatively largeQ and, in addition, has 
a string tension (the vortex mass) that depends on vortex co-ordinates through the cutoff 
in the renormalization logarithm. Then, on the large (^> £) scales we have to consider the 
vortex can not be considered as a 'relativistic' Nambu string living in the spontaneously 
broken symmetry vacuum of the dense superffuid. Given that the velocities the vortex 
reaches are on these scales always less than result, the velocity dependent part 

of the vortex self-action is dominated by that of the Magnus action. This corresponds to 
the dominance of the minimal coupling term in the momentum over the kinetic part, see 
equation ( |4T)| ) below. 

1. Co-ordinate frames 

The co-ordinates we will be using to describe the vortex self-action are defined by local 
right-handed basis vectors on the string. We can, for example, choose them to be the triad 

e» = e 2 xf, e 2 (a) = --X", f=^-X', (40) 

7 

that is, the binormal, negative normal and tangent unit vectors of the line. They are related 



by the Serret-Frenet formulas [38] to fundamental invariant properties of the line 



-^r' = -^-e 2 , -^=e 2 ' = ^ t-T?!, e x ' = Te 2 , (41) 

where R~ 1 (a) = \X"/j\ is the curvature (R c the curvature radius), and T(cx) the torsion 
of the line. For a three-dimensional superfluid a typical equilibrium line configuration is a 
circular vortex, for which e\ = ez, e 2 = en, f = e$, in a conventional cylindrical co-ordinate 
system with the normalised triad above. 

If it appears more convenient for representation purposes, we will use a more general sys- 
tem of non-normalised co-ordinate basis vectors, which are defined formally by the deriva- 
tives [0 

e a = d/dX a , e a = d/da (a =1,2), (42) 

— * 

in our non-curved Euclidean embedding space simply acting on the position vector X, 
thereby creating the basis. We will freely use the more convenient system, indicating the 
type of system by indices a = 1, 2 for the first normalised triad and a = 1 ,2 for the latter 
general co-ordinate basis. 



4 Whereby it is meant that the characteristic velocity at the core circumference, vl = k/2tt^, is 
much less than the speed of sound, if we again take £ ~ o~lj- 
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2. The self-action 

After this preparatory work of fixing conventions, we write down the vortex self-action as 
a sum of its static part and one quadratic in derivatives of perturbations from the equilibrium 
string configuration: 

S Be] i[Q{t,a)} = - J j dtday^M \^M c 2 s - \M^ 2 + ^<?' 2 } • (43) 

The displacements Q perpendicular to the line have to be small. Only in this case of small 
perturbations from an equilibrium configuration the equations of motion of interacting vortex 
segments obey the Hamiltonian structure we wish to derive below |3£j. In particular, self- 
crossings of a line bending back on itself have to be excluded. The first terms in the self 
action then represent 'nonrelativistic' terms of static and kinetic energy in the Lagrangian 
for a particle, weighted with the local differential string arc length, and integrated over the 
length of the string. 

The cutoff related elastic coefficient a (parameterized by the longitudinal as well as the 
transversal string core structure), depends on t, a in general, just as the vortex effective 
hydrodynamic mass does. In accordance with our 'nonrelativistic' treatment, we do not 
consider short wavelength perturbations, where this dependence becomes significant, and 
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take a = M c 2 = E se if. This specific choice in the elastic energy of the string is related 
to a cutoff choice in the localized self induction approximation of classical hydrodynamics 
iOfl . It corresponds to the assumption that in the long wavelength limit we are using, only 



massless sound excitations propagating along the string can survive^] It also naturally 
accounts for the fact that a must remain finite in the incompressible limit c s — > oo (whereas 
the hydrodynamic mass has to vanish). 

For a simple straight line in ^-direction, we can write for (K3T), Fourier analysing the co- 
ordinates Q = (l/(27r) 2 ) J duj J dkexp[—i(uit — kz)]Q(u,k), and neglecting the dependence 
of Mo on the co-ordinates, 

S seU /M c 2 s = - [ [dtdz + -i- / fdwdk \q 2 {u, k) { (oo/c s ) 2 - k 2 } . (44) 



,2tt) 2 J J 2 

If the wiggles of the line are occuring on scales <C R c , this should also hold for other shapes 
of line by summing up contributions of approximately straight segments. 



3. The total vortex action and momentum 

To the self-action, we have to add the interaction with the background. Using the 
conventions (^), ([31~D , 



S M = 1s jjdtda (po^ + kjX'iX^ (45) 

= q J j) dtda (-a° + OiX^j . 

The gauge potential (ipi, 6y) is defined to belong entirely to the background: 

Ofiv{x) = ttt-tj / d^kb^ityexplik^] . (46) 
\. Z7[ ) Jk<k 

The cutoff (fc At )o = (— ^/c s ,fcj) indicates the separation line between what we consider as 
being a phonon, which is integrated out in the self-action, and what we lump together into 
a time dependent and inhomogeneous background. In what follows, we will neglect the 
remaining phonon fluctuations and approximate the background to be an incompressible 
superfluid. The separation line is fixed in a quite natural way by the infrared cutoff in the 



vortex energy fl34), |fc | = expC/8i? c 



Summation of ([45]) and (|43f) yields the total vortex action 

Sy = S se li + Sm ■ (47) 



5 This is true given that we consider the kinetic and elastic parts isolated from the static energy. 
Neglecting the kinetic energy, we can obtain the Kelvin modes [p9[| . Apart from the fact that the 
vortex has a certain core structure (such that an ultraviolet cutoff is present), they are not related 
to compressibility, i.e. they exist even for infinite c s . 
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Taking the functional derivative of Sy after the vortex velocity, we arrive at the following 
expression for the vortex momentum per a-length interval of the vortex line (X = Q): 



P = p inc + p kin = qa + Mq^X . 



(48) 



It has, as expected, the same appearance as the canonical momentum for a nonrelativistically 
moving particle of charge q, which is subject to an external vector potential a. We separated 
the momentum into a part due to the (incompressible) background P mc and a kinetic (vortex 
matter) part P kin . 

Corresponding to the dominance of the static energy in the self-action, the ratio of the 
contributions to the momentum is in order of magnitude 



P 



kin 



P ; 



\x\ 



\X\c. c s 



(49) 



This depends on |X|/c a as well as k/(c s \X\) (= 0(£/\X\) in helium II). Both quantities 
are necessarily <gC 1 if the vortex is to be described within the hydrodynamic formalism we 
presented. We used here that \bij\ is of order \X k \ in a Cartesian frame (see below), and 
neglected the dependence of the self-energy logarithm on the vortex co-ordinates (which 
multiplies the right-hand side of the equation above). 



C. Gauge dependence of the vortex momentum 



The momentum (|48|) is gauge dependent through the choice for the vector potential a. 
The relation (|20| ) determines a in terms of the original field bij. This equation is, in turn, 
equivalent to the last expression in (|30"D in its nonrelativistic version, 



rot a = S = -poX' . (50) 

In a Cartesian frame, a possible solution for is b^ = — (1/3) poriijkX k (it is, for example, 
used to represent to vortex velocity dependent part of the Magnus action in j42|). This is 
an isotropic solution, which is convenient to describe the purely spatial part of the gauge 
2-form in the bulk superfluid. In the presence of the vortex, however, we choose another 
solution more appropriate to a natural basis on the string. For a singly quantized vortex, the 
momentum components P a in the basis of e a , e a = X' (a = 1 , 2), with determinant 
g = det[e*j • e}], obey (T§ 



d 2 P i r-d 1 Pr = hp ^g-, (51) 

as follows from ([50]). A factor N v in front of the right hand side enters for vortices of 
arbitrary winding number. 



An isotropic solution of ([51]) in Cartesian co-ordinates reads P mc = (l/2)hpoX x X' . 



The simplest possible solutions are obtained if we gauge one of the P- nc, s to zero: 
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pi- = ^ pine = hpQ J dx 2^ 

Pf C = => PjT = -hp, J dX 1 ^. (52) 

It is stressed that the canonical local momentum P mc is a gauge object, and not nec- 
essarily identical with a physical momentum of the vortex. A connection to a physical 
momentum can be established by using the first choice in the equation above for a circular 
vortex [e 1 = ez, eg = en, e a = X' = e$ and y/g = r), and integrating P] nc = hp Q J dX 2 yjg 
over $ = a = [0 : 2tt]. We then obtain that p Kelvin = hp 7iR 2 ez equals the Kelvin mo- 
mentum. For a superfluid the Kelvin momentum is just the surface area of the vortex 
times the bulk superfluid density times Planck's quantum of action. This result is also ob- 
tained by integrating P = J J J mpov s dV = h J J J poV9 dV = poh J f9dS, and taking the 
27r-discontinuity of 9 (cf. Figure [I]) across the surface enclosed by the vortex line. 



D. Quantization of Vortex Motion 

The canonical quantization of the massive vortex string proceeds in the usual manner. We 
take as the canonical phase space the displacement vector Q together with the momentum 
P from (|iBD . If we are to quantize the vortex motion, we have to impose the following 
canonical commutation relations, written for convenience in the basis fl40D , 

[Q», P b {a')\ = m\8{a - a') , (a, b = 1, 2) 

[Q»,QV)] = (53) 

[P a {a),P b {a')} = 

It should be observed that in our 'magnetic' field, the kinetic momentum components, 
defined in (0), do not commute, as follows from the classical Poisson bracket, respectively 
from direct calculation, using the commutators above: 

[P a ki », P b ki V)] = ihq(-p )n ab 5(a - a') , (54) 

which amounts to saying that the local velocity components of the string are not both 
determined to arbitrary accuracy. This statement is analogous to that for the velocity 
components of an electron in a magnetic field (cf. ||44|| , §110), quite as it should, according 
to the analogy explained in section |11 Bl 
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FIG. 4. The direction of the incompressible part of the vortex momentum in co-ordinate space 
depends on the choice of the gauge in (50) respectively (|5T|). On every point a of the vortex line the 
momentum can point in a different direction of the local co-ordinate plane ex, eg . This direction 
is parameterized by the angle a(a). The first choice in (|52| ) corresponds to a = 0, the second to 
tt/2. 



If we neglect the part P kin as compared to the dominant P mc in ( |48[) altogether, as indi- 
cated by (|4*D|), taking again a ring vortex as the simple archetypical example, the complete 
quantum dynamics of the (undeformed) vortex is effectively one- dimensional (has only one 
independent co-ordinate and momentum component), and given by the canonical commu- 
tator 

[Z{a), P'T (a')] = hz, hp R 2 } = ih5{a - a') 



[Z, S] = i(2%p c 



(55) 



where the second line is the version of the local commutator in the first line integrated along 
the string, so that the commutator involves the Kelvin momentum of the ring. Consequently, 
S = ttR 2 represents a surface operator of the ring plane. In the limit of P km — > 0, generally, 
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phase space and configuration space merge and become indistinguishable, the momentum 
components, then, becoming functions of the co-ordinates alone. 

The momentum space counterpart of the Figure |3] is depicted in Figure f|. We see that 
on every point of the line, we are free to choose the direction of the momentum afresh, 
according to a solution of fl5"T|) with some local gauge. 



E. The Hamiltonian and the equations of motion 



The vortex moves without dissipation, that is, there is no damping force term in the col- 
lective co-ordinate equation of motion for the vortex. Hence we can write down a conserved 
vortex energy and thus a vortex Hamiltonian. In the present case, it is useful to take the 
Coulomb gauge div a = div P mc = for the representation of the Hamiltonian: 



+ q J da r-a° c + a°^j . (56) 
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We separated Coulomb and background velocity parts of the scalar potential, i.e. wrote 



a° = a 



Q c + aP u . This gives the correct factor of 1/2 in the Coulomb interaction energy with 
other vortices: The energy of these other vortices is contained in the background part of the 
energy, whereas (l/2)ga^ is the energy solely pertaining to the vortex under consideration. 
The Hamiltonian Hy gives the total energy of this particular vortex and yields its equations 
of motion. This is true provided that the background is in the limit of infinite extension, 
such that its energy change as the vortex moves, expanding and contracting, is negligible. 
The full Hamiltonian equations of motion 

% = -l,H r , ^ = -U, (57) 
dt SQ dt sp 

in which 5/5Q,5/5P formally indicate invariant functional derivatives, are taking the form 
(setting c s = 1), 

f = -v«fc/wo - fa + <) - i (p - qS ) 2 v, (^y 



(58) 



and, of course, Q = (P — qa)/(y/yMo) . They are considerably complicated by the fact that 
we admit a dependence of the terms containing Mo, 7 on the co-ordinates of the line element. 

The two (vectorial) Hamiltonian equations of motion, being of first order in time, can be 
shown to be equivalent to the second order Lagrangian equations of motion as they follow 
from the corresponding action ([47j): 
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= q (E + (J X B) , (59) 

by making use of the identities 

dP d / '-A ■ '- 

~dt = dt (^ MoQ ) +q3+ qQ(V Q ® ^ ' 

Q x rot a = Q(V ® a) - (V ® a)Q , (60) 
^=-VgQa^ + a^ -5, Scrota. 

In these equations, the expression V ® a means a second-rank tensor with components diaj 
and a vector standing to the left or right of this tensor is contracted with the first or second 
index, respectively. 

It is instructive to write down the equation of motion in the original hydrodynamic 
variables: 

I (V7ln[- • -]Q) - 1 (-Lb* • ]Q') + (i - W H ■ ■]) 

+iq' 2 v 5 (-L '"[■•■]) = (4Vr,)A"x ((j-ff.) . (6i) 

The logarithm of the self-energy ( R^ ) is abbreviated ln[- • •]. The bulk superfluid density po 
disappears in this writing. 



IV. SUMMARY 

In this paper, we have developed a general hydrodynamic formalism to describe the 
zero temperature motion of vortices in compressible, conventional superfluids, and the fields 
emanating from and interacting with them during this motion. The assumptions pertaining 
to this formalism are as follows. There exists an entity vortex, whose center of topological 
stability, the singularity of the vortex center, can be described by the spacetime embedding 
of the vortex world sheet. The intrinsic property of the vortex is the circulation, defined 
as the line integral of the momentum p, and quantized in units of Planck's quantum of 
action h by the homotopy group index iV^. In the underlying spacetime, there exists a 
conserved particle current density j(x), whose dual is the field H acting on the vortex. 
The fundamental degrees of freedom in this treatment are thus the position of the vortex 
X^ir, a) in spacetime, representing the order parameter defect location, and the current 
density j(x). In a perfect fluid, these two degrees of freedom are connected physically in the 
Magnus force action. This action is given by the linear coupling of the vorticity u> = *d A p 
and a gauge potential b, whose exterior derivative gives the field strength H = d A b. 

In the Galilei invariant, nonrelativistic fluid we are then dealing with, we have further 
separated the superfluid phase into the nonrelativistic part of the background fluid and a part 
related to small density and phase oscillations, namely sound waves. This latter part behaves 



22 



'relativistic' in that the equations governing its behaviour have pseudo-Lorentz invariance, 
that is, they are Lorentz invariant under the replacement c — > c s . Self-interaction of the 
vortex line with these 'relativistic' phonon excitations gives rise to the hydrodynamic vortex 
mass. In a dense superffuid, the 'nonrelativistic' motion of the vortex is of relevance, because 
large order parameter modulus (quantum) variations take over, long before \dX/dt\ < c s is 
attained. The full canonical structure of vortex motion obtains by adding elastic energy in 
the self-action of the vortex string, related to the localized self induction approximation. 

In the final description of the vortex as a massive, elastic string object, we have thus 
incorporated parts of the fundamental interaction of the bare vortex with the general flow 
field. The phonon fluctuations in this flow field give mass renormalization, whereas the 
interaction with adjacent line segments gives elasticity. The background flow potential which 
remains is well approximated to be coming from an incompressible fluid. This procedure 
yields equations of motion in which the role of the bare, small scale many-body quantum 
dynamics of the vortex is reduced to give an ultraviolet cutoff parameter in the mass and 
elasticity coefficients. 

The local self-action of the laboratory vortex, considering only the influence of neigh- 
bouring segments on a particular line element, has to be modified in a nonlocal manner if 
more distant parts of the line are closely approaching each other, or when nontrivial topolo- 
gies of the vortex line, knotted structures, are under consideration Given that possible 



vortex configurations are regular and topologically trivial in this sense, the equations of 



motion in subsection 1!IE can, for example, be used to calculate the Euclidean action of 



vortex tunnelling motions under the influence of external flow fields [|15| , |16 
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